Genome-Wide Association and Expression Analysis Revealed the Candidate Variants and Molecular Underpinnings of Cold-Stress Response in Large Yellow Croaker.

Large yellow croaker (Larimichthys crocea) is one of the most economically important fish in China. Recently, global climate change has caused more and more intense and extreme low temperature weathers, resulting in huge losses to the large yellow croaker industry. Therefore, it is essential to understand the mechanisms of low-temperature tolerance in large yellow croaker. Here, we conducted an integrative analysis of genome-wide association study (GWAS) and transcriptome analysis to identify candidate variants and reveal the molecular underpinning of cold-stress response in large yellow croaker. A total of 8 significant single nucleotide polymorphisms (SNPs) loci on 6 chromosomes were identified in the GWAS analysis, and 5764 (gill) and 3588 (liver) differentially expressed genes (DEGs) were detected in cold-stressed large yellow croaker, respectively. Further comparative and functional analysis of the candidate genes and DEGs highlighted the importance of pathways/genes related to immune response, cellular stress response, lipid transport, and metabolism in the cold-stress response of large yellow croaker. Our results provide insights into the cold tolerance of large yellow croaker and contribute to genomic-based selection for low-temperature-resistant large yellow croaker.

Genome-Wide Association Study Identifies Genomic Loci of Sex Determination and Gonadosomatic Index Traits in Large Yellow Croaker (Larimichthys crocea).

Larimichthys crocea is one of the traditional marine culture fishes in China, widely distributed in South China Sea, East Sea, and southern Yellow Sea. Sex dimorphism is evident in this species that females present a substantial growth strength than males, suggesting breeding females could obtain more economic benefits in L. crocea aquaculture industry. With the continuous expansion of aquaculture industry, both identifying sex-associated genome region and understanding the genetic basis underlying gonad differentiation and development matter to not only sex control aquaculture but also breeding industry. Thus, genome-wide association analysis (GWAS) of sex determination was conducted with a random breeding population of 905 individuals (including 463 females and 442 males) by ddRAD sequencing. For sex determination, 21 significant single nucleotide polymorphisms (SNPs) in chromosome (Chr) 22 were identified. Surrounding these SNPs, we founded 14 candidate genes, including dmrt1, dmrt3, and piwil2, fam102a, and odf2. The sex-associated region was narrowed down further to 2.4 Mb on Chr22 through Fst scanning and insertion-deletion (InDel) analysis. Besides, 3 SNPs in the supposed sex-determining region on Chr22 were identified as highly associated with gonad differentiation through GWAS on gonadosomatic index (GSI) in 350 males and 231 females. Because of the significant difference of GSI between females and males of L. crocea, GWAS on GSI of different genders was also conducted independently. Finally, we identified a SNP in Chr18 showing genome-wide significant association with male GSI (MGSI) and three genes axl, cyp2a10, and cyp2g1 involved in the gonadal development regulation process of aromatase. Overall, this study explored the genetic basis of sex determination mechanism and provided novel insights into gonad differentiation and development, offering solid genetic support for sex control breeding, marker-assisted selection, and marine resources conservation.

Genome-Wide Association Study of Body Shape-Related Traits in Large Yellow Croaker (Larimichthys crocea).

Large yellow croaker (Larimichthys crocea) is one of the most important cultured marine fish on the southeast coast of China. Its body shape is important for the aquaculture industry since it affects the behavior such as swimming, ingesting, and evading, as well as customer preference. Due to the greater consumer demand of small head, slender body large yellow croaker, selecting and breeding of slender individuals with the assistance of genetic markers will benefit the industry quickly. In this study, several traits were employed to represent body shape, including body depth/body length (BD/BL), body thickness/body length (BT/BL), caudal peduncle depth/caudal peduncle length (CPDLR), tail length/body length (TL/BL), and body area/head area (BA/HA). Genome-wide association study was conducted with a panmictic population of 280 individuals to identify SNP and genes potentially associated with body shape. A set of 20 SNPs on 12 chromosomes were identified to be significantly associated with body shape-related traits. Besides, 5 SNPs were identified to be suggestive associated with CPDLR and BT/BL. Surrounding these SNPs, we found some body shape-related candidate genes, including fabp1, acrv1, bcor, mstn, bambi, and neo1, which involved in lipid metabolism, TGF-ß signaling, and BMP pathway and other important regulatory pathways. These results will be useful for the understanding of the genetic basis of body shape formation and helpful for body shape controlling of large yellow croaker by using marker-assisted selection.